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OO ; 

O Supernova remnants(SNRs) are believed to be the primary lo- 

^ | cation of the acceleration of Galactic cosmic rays, via diffusive 

shock (Fermi) acceleration. Despite considerable theoretical work 
the precise details are still unknown, in part because of the diffi- 
^ | culty in directly observing nucleons that are accelerated to TeV 

energies in, and affect the structure of, the SNR shocks. How- 
ever, for the last ten years, X-ray observatories ASCA, and more 
recently Chandra, XMM-Newton, and Suzaku have made it possible 
' to image the synchrotron emission at keV energies produced by 

^ ■ cosmic-ray electrons accelerated in the SNR shocks. In this article, 

we describe a spatially-resolved spectroscopic analysis of Chandra 
observations of the Galactic SNR Cassiopeia A to map the cutoff 
frequencies of electrons accelerated in the forward shock. We set 
upper limits on the diffusion coefficient and find locations where 
particles appear to be accelerated nearly as fast as theoretically 
possible (the Bohm limit). 

Supernova remnants (SNRs) have been established as the leading can- 
didate for the acceleration of cosmic rays.pP|2] It has been shown that the 
mechanism of diffusive shock acceleration in SNR shocks coupled with some 
understanding of Galactic transport effects can in theory produce the ob- 
served power-law spectrum of cosmic rays. jSJ-^Ol- This model works at least 
up to the "knee" of the cosmic ray spectrum near 5 x 10 15 eV, and possibly 
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all the way to the "ankle" near 3 x 10 18 eV. jllj The charged particles scatter 
off of irregularities in the magnetic field, increasing their momentum by a 
fraction of the shock velocity v/c with each round-trip shock crossing, jllj 
Theoretical work in the last several years has suggested that the process is 
significantly nonlinear . 1 1 Uj - [To] Higher energy particles have larger diffusion 
lengths and therefore sample a greater change in velocity and compression 
ratio across the shock than lower energy particles. This effect introduces a 
decrease in the particle distribution's spectral index and a flattening of the 
spectrum with increasing energy. [Til] A smoothing effect on the structure of 
the shock results, predominantly caused by the ions; electrons are almost 
"test particles." At high energies the electron spectrum is expected to have 
the same (curved) spectral shape as the proton spectrum. [TU] 

This correspondence of electron and proton spectra is significant because 
it is extremely difficult to directly observe the acceleration of cosmic-ray 
protons in SNR shocks. Their presence is inferred indirectly by detecting 
7T° decay 7-rays produced from collisions of the cosmic rays with gas par- 
ticles in the ambient medium. Atmospheric Cerenkov telescopes have ob- 
served gamma rays from some supernovae: Cas A with HEGRA 16*] |17j. 
RX J1713.7-3946 (G347.3-0.5) with CANGAROOQ"""| d and HESSjU], 
RX J0852.0-4622 with HESS"2"Q and CANGAROO-II"""""", and G 0.9+0.1 
with HESS 23. Current technology limits the spatial resolution of these in- 
struments to arcminutes and remnants the size of Cas A may appear point- 
like. In addition, TeV gamma rays can also be produced by inverse-Compton 
scattering of cosmic microwave background photons with the accelerated 
cosmic-ray electrons in SNRs. Consequently, it is difficult to precisely and 
unambiguously image proton and ion acceleration in the SNR. In contrast, 
the electrons which are accelerated to TeV energies produce synchrotron 
radiation at keV energies which is relatively easily imaged at arcsecond res- 
olution. The electrons are not expected to be accelerated all the way to 
"knee" energies, as at some point the synchrotron losses equal the accelera- 
tive gains on shock crossing and the spectrum cuts off. However, it is still 
possible to set a lower limit on the efficiency of the acceleration. 

An analysis of ASCA observations of SN 1006 provided the first evidence 
of X-ray synchrotron emission in an SNR. |24| . The approximately arcsecond 
spatial resolution over much of the energy band (0.5 - 10 keV) of the Chan- 
dra High Resolution Mirror Assembly-the best of any X-ray telescope which 
has flown j25] -has enabled detailed studies. For example, Lazendic et al. 26 
studied Chandra observations of G347.3-0.5 to localize cosmic ray accelera- 
tion in filaments. Since the X-ray spectra of the youngest Galactic supernova 
remnants (Cas A, Tycho, and Kepler) are dominated by thermal emission, 
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the ability to spatially resolve the nonthermal contribution is especially im- 
portant. Using Chandra observations, Warren et al. [27] demonstrated that 
acceleration in Tycho occurs at the forward shock. 

Chandra's Advanced CCD Imaging Spectrometer's (ACIS) spectral res- 
olution (E/AE ~ 10 — 60) allows extraction of the spectrum of the emission. 
The combined capabilities of the mirror and CCDs allows us to determine 
if a nonthermal spectrum is present and where it is located. We: 

1) developed a method to precisely identify the regions primarily radi- 
ating nonthermal X-ray synchrotron emission in a remnant dominated by 
thermal X-ray radiation from supernova ejecta; 

2) fit the spectra of thousands of individual regions dominated by non- 
thermal emission in Cas A with a realistic synchrotron model, to map the 
synchrotron cut-off frequency around the outer shock; 

3) determined the ratio of the estimated upper limit on the diffusion co- 
efficient, k, in each region to the Bohm coefficient, Kb, and thereby mapped 
the lower limit of the efficiency of acceleration in a SNR. In several areas, 
diffusion appears to be occurring near the Bohm limit. The maximum effi- 
ciency occurs in the northeast outer shock, just to the side of the jet, where 
k/kb < 2.1. 

Mapping the spectral properties of the Cas A SNR 

Exploiting the full capabilities of Chandra's spatial resolution, ACIS's spec- 
tral resolution, the scripting functionality of the Interactive Spectral Inter- 
pretation System (ISIS)[2H], and distributed computation, we produced high 
quality, high resolution maps of the spectral parameters of the remnant. We 
used direct spectral fitting, not energy selection methods. This process al- 
lowed us to map parameters such as temperature and interstellar absorption 
column, which is not possible with other methods. 

We used the nine individual pointings from the million second obser- 
vation of Cas Aj22j and two additional archival observations to take the 
deepest look to date at a SNR with Chandra. Figure ^ shows a combined 
counts image color-coded by energy. To perform the spatially-resolved spec- 
tral analysis, we divided the sky around Cas A into a square grid of sky- 
coordinate region centers separated by approximately 1" (2 ACIS pixels). 
Cas A was covered by 104,393 regions. Each region was adaptively sized to 
contain a minimum of 10,000 counts. The region sizes ranged from 1" x 1" 
to 7" x 7" on source. Regions were allowed to overlap; therefore, fits to 
regions larger than 1" x 1" are not independent. 

The position and roll angle of Cas A on the CCD varies in the eleven 
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pointings. We applied standard processing and filtering to the data for each 
pointing and then some custom techniques. [3*0] We re-projected the data 
to a common sky-coordinate plane, but maintained separate files for each 
pointing. To speed the later extraction of spectra, we histogrammed the 
spectral data and determined the mean chip coordinates and the exposure 
time for each sky location. This allowed us to separate the slowly varying 
effect of the effective area, which we calculated on a 32 pixel by 32 pixel 
grid, and the quickly varying effect of bad pixels and the exposure time, 
which we calculated on a pixel to pixel basis. Since our modified effective 
area depends only on chip location, we used one library of the effective area 
function for all of the pointings corresponding to a given observation epoch. 
Similarly, we pre-computed a library of energy response matrices on the 
same 32 pixel by 32 pixel grid. 

Using ISIS, the same spectral model was fit jointly to the set of eleven 
spectra for each region. Each spectrum was associated with its own ef- 
fective exposure, response functions, and optionally, background spectrum 
extracted from the same pointing. The background spectrum was not sub- 
tracted from the source spectrum, but was instead added to the model. The 
background is dominated by photons from Cas A in the wings of the point 
spread function of the telescope and in the CCD readout streak. For com- 
putational simplicity, the background for every source spectrum was drawn 
from the same representative sky region. 

The analyses were performed in two phases. In Phase 1, we used a sim- 
ple fitting function of fifteen Gaussian lines for specific O, Ne, Mg, Si, S, Ar, 
Ca, Ti and Fe line complexes (representing He-like K Q for all, He- like 
for Si and S, H-like K Q for Si, and L lines for multiple Fe ions); an interstel- 
lar absorption column; and a bremsstrahlung continuum. This particular 
model was fit without background spectra (because it fits the background) 
to the 104,393 regions covering the entire remnant. The joint best-fit spec- 
tral parameters were recorded for each region and used to create FITS im- 
ages of the spectral properties. Figure El shows the fit temperatures of the 
bremsstrahlung continuum. Structures concentrated at the outer edges and 
in part of the center show a significantly higher temperature than the rest of 
the remnant. As shown in Figures [2 and 01 these regions show little evidence 
of line emission and are correlated with anomalously high absorption column 
values and with the high energy (4-6 keV) continuum-dominated filaments 
in Figure ^ There appears to be some genuine elevation of the absorption 
column in the western region of the remnant, where there is evidence that 
the remnant is obscured by a molecular cloud 31 . At the outermost edges, 
the filaments have been associated with the outer shock. (223- Collectively, 
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these clues suggest that the model used for the first phase of the analysis 
is inappropriate in the anomalously high temperature regions because the 
emission is nonthermal. While some of the high energy X-ray emission (10- 
60 keV) |33| from Cas A may be produced by nonthermal bremsstrahhmg 
|34|-|36|. there is general agreement that the X-ray emission observed from 
forward-shock filaments is produced by synchrotron emission. Therefore, in 
Phase 2, we fit the spectra of regions identified by this method with a model 
composed of realistic synchrotron emission, an absorption column, and a 
background spectrum. 

Creating a cut-off frequency map 

Having found a means to identify the regions dominated by nonthermal 
emission in a remnant whose X-ray emission is globally dominated by ther- 
mal emission, we studied the properties of particle acceleration. We used 
the kT results to pick 10,857 regions-the top 10% of fitted temperatures 
(kT > 2.6 keV). This criterion clearly selected areas centered on and fully 
containing the outer filaments but not excessive surrounding areas or iso- 
lated peaks. Since the fit cut-off frequencies span a range of two orders 
of magnitude, we believe we found all of the regions of interest without 
introducing selection effects in the distribution of detected cut-offs. 

Our group developed a synchrotron spectrum model in which the effects 
of spectral "curvature" from the back-pressure of the accelerated cosmic rays 
are modelled using an electron spectrum whose spectral index is a linear 
function of the logarithm of the momentum. |37| . The core of the model is 
the electron distribution function, 



in which n is the electron number density, c is the speed of light, p = ^ymv 
with m the electron mass and v the velocity of a particle, A e is the number 
density at p =1 GeV/c, T is the differential spectral index at p =1 GeV/c, e 
is the exponential cut-off (maximum) energy, and a is the spectral "curva- 
ture." The parameters of the synchrotron model are the normalization N s 
(a function of A e , distance, and the emitting volume), the total magnetic 
field -Btoti r> a i and e. 

The cut-off energy e, the normalization of the spectrum and the value 
of an interstellar absorption column were allowed to be free parameters. 
We froze the curvature parameter to a value of 0.06 based on fits to radio 
and infrared synchrotron data for selected regions of Cas A [HE]- A value 




dpc 



GeV 




5 



of a = 0.05 was obtained by our group in fits to SNf006 39 . As there 
are no published values for the radio spectral index of Cas A at the outer 
shocks, which are at the edge of a diffuse plateau of radio emission, we set 
r = 2.54, using the average radio spectral index for the entire remnant. _40 
Comparison with the index in the plateau region just inside the outer X-ray 
filaments from unpublished radio spectral index maps suggests this value is 
reasonable. Small variations in the value of T did not significantly affect our 
final results. Changing the value of V by ±0.1 changed the cut-off frequencies 
to values within our confidence limits and preserved overall trends. 

Some assume the widths of the filaments are limited by synchrotron 
losses and estimate the strength of the magnetic field in the filaments of 
Cas A to be « 100 ^GjHn] to as 500 //G[3I]. Others argue the sizes of the 
filaments in SNRs are determined by the widths of regions of enhanced mag- 
netic fields at the forward shock, limited by strong turbulent damping of the 
field to physical sizes of order 10 16 - 10 17 cm|42j. This range corresponds to 
0.2 to 2" at the distance of Cas A, similar in size to the observed width. The 
maximum frequency of emission from electrons is related to the maximum 
energy and the magnetic field as 

It may not be possible to uniquely determine e and B tot based on solely 
on analyses of the X-ray spectra. Therefore, we froze -Btot = 1000/iG (the 
"equipartition value" for Cas A[HB])- As a result, the cut-off energy itself 
is not well determined, but the cut-off frequency is well determined by the 
X-ray data. 

We present our map of v c in Figure Cut-off frequencies in the forward 
shock vary at least from 5 xlO 16 Hz to 9xl0 17 Hz. The 90% confidence 
intervals have been computed for all of the cut-off frequencies. The average 
uncertainties range from about ± 20% at 5 x 10 16 Hz to ± 40% at 9x 10 17 Hz. 
Therefore, the cut-off energy and/or the magnetic field vary considerably 
from one region to another. Variations in the cut-off frequency have been 
reported for a few other remnants 26 j . [13]- 145| . The spectrum of a 2" x 2" 
region along the shock in the northeast is shown in Figure El This region 
has the highest cut-off frequency (9xl0 17 Hz). 

Constraining the electron diffusion coefficient 

As recently reviewed by VinkjjUi], some analyses are based on the assump- 
tion that particle acceleration is efficient: the particle diffusion coefficient is 
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as small as the Bohm limit |47| |48j. Here, we use the spectral results to con- 
strain the efficiency of acceleration of the highest energy electrons at each 
location dominated by synchrotron radiation. 

Assuming that electrons at the cut-off energy cannot experience syn- 
chrotron losses in excess of accelerative gains in one cycle through the shock, 
one can derive an upper limit on the ratio of the electron diffusion coefficient 
to the Bohm coefficient. |39|. The ratio is a function of the cut-off frequency, 
the shock velocity u\ and the compression ratio r. At E = e, Kb = e/3eB 
in MKSA and 



k_ < 9e mc fuf 
kb ~ 2e 2 v c sin 6 



0.15 J V5000 km s' 1 J \9 x 10 17 Hz, 

Here eo is the permittivity of free space, and / = (r — l)/((r)(r + 1)). 

To constrain the electron acceleration efficiency, we used a value of 
5000 km/s for the shock velocity of Cas Aj^H] ; and assumed the case of 
a strong, unmodified shock (/ = 0.15, r = 4). For the highest cut-off fre- 
quencies on the outer shock of Cas A, we obtained k/kb = 2.lio.7- The 
variations on the outer rim range from 0.5 to 9xl0 17 Hz, corresponding 
to k/kb ~ 36 to 2. Our upper limits are smallest in certain parts of the 
north, northeast, and southeast shocks, strongly suggesting efficient particle 
acceleration occurs in at least these regions. 



Conclusions 



We developed a method to analyze the emission from extended X-ray sources. 
We combined multiple data sets retaining pointing-specific calibrations and 
backgrounds, automated spectral extraction, and fully exploited the spa- 
tial and spectral capabilities of Chandra. Analyzing the emission from 
the young Cassiopeia A supernova remnant, we discovered the apparent 
bremsstrahlung temperature (kT) of regions can be used to identify individ- 
ual regions which are dominated by nonthermal emission in a SNR globally 
dominated by thermal emission. For Cas A, we found that the synchrotron 
emission is clearly correlated with the forward shock, as well as some areas 
in the center, which may or may not be the forward shock seen in projection. 
Fitting the nonthermally dominated regions with a synchrotron model, we 
mapped the cut-off frequency of the synchrotron spectrum associated with 
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the highest energy electrons accelerated in the forward shock. From the 
cut-off frequency and shock velocity, we constrained the efficiency of the 
acceleration of cosmic-ray electrons. 

We present the first map of the upper limits on the particle diffusion 
coefficients in a SNR. Our results suggest that at least in several locations, 
acceleration occurs nearly as fast as theoretically possible (the Bohm limit). 
In fact, electrons could also be diffusing in the Bohm limit in other parts of 
the shock if the shock velocity is lower; if the compression ratio is greater 
than 4 — an enhancement quite possible in a modified shock; or if a process 
other than synchrotron radiation limits the maximum energy. TeV 7-rays 
have been detected from Cas A^HICZI- Modeling suggests these could be 
formed either by the prompt decay of 7r°'s created in collisions of cosmic-ray 
protons accelerated in the SNR with ambient gas; or, by inverse-Compton 
scattering of cosmic microwave background photons off of the accelerated 
electrons. Although TeV telescopes cannot yet localize the 7-ray source po- 
sition with arcsecond accuracy, or determine unambiguously the fraction of 
7T° and inverse-Compton produced 7-rays from the spectrum, the expecta- 
tion that cosmic-ray electrons and protons are accelerated together means 
that electron synchrotron emission is an important tracer of the particle ac- 
celeration. Detailed maps of the X-ray synchrotron emission, such as those 
we created with Chandra, are rapidly improving our understanding of the 
big picture of cosmic-ray acceleration in SNRs. 
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Figure 1: The Cassiopeia A supernova remnant as observed by Chandra. 
Red: 0.5 - 1.5 keV: O-K, Ne-K, Mg-K, and Fe-L emission. Green: 1.5-2.5 
keV: Si-K and S-K. Blue: 4.0-6.0 keV: high energy continuum (some energy 
ranges have been excluded, including Fe-K). This energy-colored photon 
counts image was created from 1.08 Ms of ACIS data. The colors are indi- 
vidually log-scaled to bring out the fine structure in the remnant. The bulk 
of the X-ray photons are thermal radiation from the ejecta which has passed 
through the hot reverse shock. Note in particular the wispy blue filamen- 
tary structure of the forward shock, dominated by high-energy synchrotron 
emission from electrons accelerated at the forward shock. 
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Figure 2: Cas A bremsstrahlung temperature map. The image is of the 
fit bremsstrahlung continuum from the Phase 1 analysis of the remnant. 
We fit the spectrum of 104,393 regions with a model that includes 15 Gaus- 
sian lines, interstellar absorption, and a thermal bremsstrahlung continuum. 
Note the correlation between the high-temperature features in this image 
and the blue, high-energy continuum-dominated filaments in Figure ^ The 
anomalously high temperatures in this image indicate nonthermal emission 
in the SNR, which is likely synchrotron radiation from electrons accelerated 
to TeV energies. 
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Figure 3: Comparison of a continuum-dominated spectrum (blue) from a 
2" by 2" region on the northeast shock (a=23:23:42.2, 5=58:50:26.4) shows 
it has a different spectral shape than spectra from same-sized regions dom- 
inated by line emission (Si, S, Fe, O, Ne, Mg, Ar, Ca). The green spectrum 
is taken from a region in the northern lobe (a=23:23:52.1, (5=58:50:33.3) 
featuring strong K-line emission from silicon and sulfur. The red spec- 
trum is from the fingerlike projections in the southeast lobe (a=23:23:44.0, 
5=58:48:04.7) dominated by iron L and K emission. The exposure time is 
1.08 Ms, and for display the spectra are binned to contain at least 10 counts 
per bin. 
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Figure 4: Synchrotron cut-off frequencies and upper limits on the electron 
diffusion coefficient in the Cassiopeia A supernova remnant. A synchrotron 
model which includes curvature in the electron distribution spectrum was 
used to fit the spectra of 10,857 regions identified from the kT map as likely 
sites of nonthermal emission (see Figure EJ) because they have anomalously 
high bremsstrahlung temperatures (i.e. kT > 2.6 keV, Figure [2J- The fitted 
cut-off frequencies and estimates of the shock velocity and compression ratio 
are used to determine an upper limit on the ratio of the electron diffusion 
coefficient (R) to the Bohm diffusion coefficient (kb)- The highest fit cutoff 
value occurs at the 2" x 2" region indicated by the green cross and arrow; 
the fit is shown in Figure El A ratio close to 1, as seen in the yellow and white 
regions of the map, indicates that acceleration must be occurring almost as 
fast as theoretically possible (in the Bohm limit). 
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Figure 5: The nonthermal continuum spectrum of the 2" by 2" region on 
the northeast shock (a=23:23:42.2, 5 = 58:50:26.4) which shows the highest 
cut-off frequency, 9 xlO 17 Hz. The blue histogram is the combined data; 
the red histogram is the model plus combined background. The apparent 
(thermal) structure between 1 and 3 keV is primarily background. Only 
every third error bar has been plotted to improve clarity. The exposure 
time is 1.08 Msec and for fitting the spectrum has been binned to contain a 
minimum of 100 counts per bin. 
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